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Abstract. I review recent progress in finite temperature lattice calculations, 
including the determination of the transition temperature, equation of state, screening 
of static quarks and meson spectral functions. 
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O ■ 1. Introduction 

Ch . 

One expects that at sufficiently high temperatures and densities the strongly interacting 
matter undergoes a transition to a new state, where quarks and gluons are no longer 
confined in hadrons, and which is therefore often referred to as a deconfined phase or 
\^ ', Quark Gluon Plasma (QGP). The main goal of heavy ion experiments is to create such 
form of matter and study its properties. We would like to know at which temperature 
the transition takes place and what is the nature of the transition as well the properties 
of the deconfined phase, equation of state, static screening lengths, transport properties 
i etc. Lattice QCD can provide first principle calculation of the transition temperature, 
equation of state and static screening lengths (see Ref. [H 12]) for recent reviews. 
^ I Calculation of transport coefficients remains an open challenge for lattice QCD (see 
discussion in Ref s. OH]). 

One of the most interesting question for the lattice is the question about the nature 
of the finite temperature transition and the value of the temperature T^. where it takes 
place. For very heavy quarks we have a 1st order deconfining transition. In the case of 
QCD with three degenerate fiavors of quarks we expect a 1st order chiral transition for 
sufficiently small quark masses. In other cases there is no true phase transition but just 
a rapid crossover. Lattice simulations of 3 fiavor QCD with improved staggered quarks 
(p4) using A^^ = 4 lattices indicate that the transition is first order only for very small 
quark masses, corresponding to pseudo-scalar meson masses of about 60 MeV [lO]. A 
recent study of the transition using effective models of QCD resulted in a similar estimate 
for the boundary in the quark mass plane, where the transition is 1st order |8J. This 
makes it unlikely that for the interesting case of one heavier strange quark and two 
light u, d quarks, corresponding to 140 MeV pion, the transition is 1st order. However, 
calculations with unimproved staggered quarks suggest that the transition is 1st order for 
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pseudo-scalar meson mass of about 300 MeV [H]. Thus the effect of the improvement 
is significant and we may expect that the improvement of flavor symmetry, which is 
broken in the staggered formulation, is very important. But even when using improved 
staggered fermions it is necessary to do the calculations at several lattice spacings in 
order to establish the continuum limit. Recently, extensive calculations have been done 
to clarify the nature of the transition in the 2+1 flavor QCD for physical quark masses 
using Nt = 4, 6, 8 and 10 lattices. These calculations were done using the so-called 
stout improved staggered fermion formulations which improves the flavor symmetry of 
staggered fermions but not the rotational symmetry. The result of this study was that the 
transition is not a true phase transition but only a rapid crossover [T2| . New calculations 
with stout action indicate that only for quark masses about ten times smaller than the 
physical quark mass the transition could be first order pL3j. Even-though there is no 
true phase transition in QCD thermodynamic observables change rapidly in a small 
temperature interval and the value of the transition temperature plays an important 
role. The flavor and quark mass dependence of many thermodynamic quantities is 
largely determined by the flavor and quark mass dependence of Tg. For example, the 
pressure normalized by its ideal gas value for pure gauge theory, 2 flavor, 2+1 flavor 
and 3 flavor QCD shows almost universal behavior as function of T/T^ [9]. 

The chiral condensate (ipip) ^ind the expectation value of the Polyakov loop (L) 
are order parameters in the limit of vanishing and infinite quark masses respectively. 
However, also for finite values of the quark masses they show a rapid change in vicinity 
of the transition temperature. Therefore they can be used to locate the transition 
temperature. The fluctuations of the chiral condensate and Polyakov loop have a peak 
at the transition temperature. The location of this peak has been used to define the 
transition temperature in the calculations with p4 action on lattices with temporal 
extent A^,- = 4 and 6 for several values of the qurk mass [7j. The combined continuum 
and chiral extrapolation then gives the value Tc = 192(7) (4)MeV. In this calculations 
the lattice spacing has been fixed by the Sommer parameter tq = 0.469(7)fm |16j . 
The last error in the above value of Tc corresponds to the estimated systematic error 
in the extrapolation. Recently the transition temperature has been determined using 
the so-called stout staggered action and Nt = 4, 6, 8 and 10 [19]. The deconfinement 
temperature has been found to be 176(3) (4) MeV |I9]. The central value is considerably 
smaller than the one obtained with p4 action but taking into account the errors the 
deviation is not very significant. The authors of Ref. [19] use a different definition 
of the chiral susceptibility which resulted in the chiral transition temperature of 
Tchirai = 151(3)(3)MeV. Using the chiral susceptibility defined above would result in a 
larger value of the transition temperature. Although the continuum extrapolation using 
only Nt- = 4 and 6 is not completely reliable new calculations by HotQCD Collaboration 
using Nr = 8 lattices give quite similar results for the Polyakov loop, chiral condensate 
and strangeness susceptibility as earlier calculations with N^- = 4 and 6 lattices both for 
the p4 and asqtad action [20] . 

Lattice calculations of equation of state were started some twenty years ago. In the 
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Figure 1. The interaction measure calculated (left) and the entropy density (right) 
for the p4 action 



case of QCD without dynamical quarks the problem has been solved, i.e. the equation 
of state has been calculated in the continuum limit [21]. At temperatures of about 4Tc 
the deviation from the ideal gas value is only about 15% suggesting that quark gluon 
plasma at this temperate is weakly interacting. Perturbative expansion of the pressure, 
however, showed very poor convergence at this temperature [22] . Only through the use 
of new re-summed perturbative techniques it was possible to get agreement with the 
lattice data [23l IMl EH]. To get a reliable calculation of the equation of state on the 
lattice, improved actions have to be used [261 [27]. Recently equation of state have been 
calculated using p4 and asqtad improved staggered fermion actions [28l [29]. In lattice 
calculation the basic thermodynamic quantity is the trace of the energy momentum 
tensor. , often refered to as the interaction measure e — 3p. This is because it can be 
expressed in terms of expectation values of gauge action and quark condensates (see 
discussion in Ref. [29]). All other thermodynamic quantities, pressure, energy density 
and entropy density s = {e + p) can be obtained from it using integration 

p(T)_p(7o)^ r ^^, e{T')-3p{r) 

^4 7^4 J^^^ ^ ) 

The value of Tq is chosen to be sufficiently small so that it corresponds to vanishing 
pressure to a fairly good approximation. In Fig. [1] I show the interaction measure 
form the new calculations with p4 action on = 4 and 6 lattices [29j. At highest 
temperatures calculations with A^.^ = 8 lattices have also been performed. As expected 
because of use of the improved action the difference between the A^.^ = 4 and A^ = 6 
results is small. In this figure I also show the entropy density which raises rapidly in the 
tempearture region 180 — 200 MeV. At high temperature it is only 10% below the ideal 
gas limit in agreement with expectations from improved perturbative calculations. The 
results from calculations with the asqtad action [28] are also shown. These calculations 
agree very well with the p4 results providing further evidence that the cutoff effects are 
small. 
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2. Spatial correlation functions 

To get further insight into properties of the quark gluon plasma one can study different 
spatial correlation functions. One of the most prominent feature of the quark gluon 
plasma is the presence of chromoelectric (Debye) screening. The easiest way to study 
chromoelectric screening is to calculate the singlet free energy of static quark anti-quark 
pair (for recent review on this see Ref. [30]), which is expressed in term of correlation 
function of temporal Wilson lines 

exp(-Fi (r,T)/T) = TT{W{r)W\0)). (2) 

L = TtW is the Polyakov loop. This quantity is also useful to study quarkonia 
binding at high temperatures [3ll |32l EHl |3ll |35l EH]- In purely gluonic theory the 
free energy grows linearly with the separation between the heavy quark and anti-quark 
in the confined phase. In presence of dynamical quarks the free energy is saturated 
at some finite value at distances of about 1 fm due to string breaking [301 EH 138] . 
Above the deconfinement temperature the singlet free energy is exponentially screened, 
at sufficiently large distances [38], i.e. 

Fi(r,T) = F^iT) - l^exp(-m^(T)r). (3) 

The inverse screening length or equivalently the Debye screening mass niD is 
proportional to the temperature. In leading order of perturbation theory it is = 
\/l + Nf /3g{T)T. Beyond leading order it is sensitive to the non-perturbative dynamics 
of the static chromomagnetic fields. The Debye screening mass has been calculated in 
pure gauge theory {Nj = 0) [38] and in 2 fiavor QCD {Nj = 2) [39j as well as in 2+1 
fiavor QCD [40j. The temperature dependence of the lattice data have been fitted with 
the simple Ansatz motivated by the leading order result : m£)(T) = A^yl + Nf /3g{T)T. 
Here g{T) is the two loop running coupling constant. This simple form can fit the data 
quite well if A ~ 1.4 — 1.6 ^40j. Thus the temperature dependence as well as the 
fiavor dependence of the Debye mass is given by perturbation theory. We also see that 
non-perturbative effects due to static magnetic ffelds significantly effect the electric 
screening, resulting in about 40% corrections. However, the non-perturbative correction 
is the same in full QCD and pure gauge theory. 



3. Spectral functions 

Information on hadron properties at finite temperature as well as transport coefiicients 
are encoded in different spectral functions. In particular the fate of different quarkonium 
states in the quark gluon plasma can studied by calculating the corresponding 
quarkonium spectral functions. On the lattice we can calculate correlation function 
in Euclidean time. The later is related to the spectral function via integral relation 
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Given the data on the Euchdean meson correlator G{t, T) the meson spectral function 
can be calculated using the Maximum Entropy Method (MEM) [51]. For charmonium 
this was done by using correlators calculated on isotropic lattices [l2l |13] as well as 
anisotropic lattices [HI US], |16] in quenched approximation. It has been found that 
quarkonium correlation function in Euclidean time show only very small temperature 
dependence [^ B6] . In other channels, namely the vector, scalar and axial- vector 
channel stronger temperature dependence was found [l3l |16], especially in the scalar 
and axial-vector channels. The spectral functions in the pseudo-scalar and vector 
channels reconstructed from MEM show peak structures which may be interpreted as a 
ground state peak [HI US, US] . Together with the weak temperature dependence of the 
correlation functions this was taken as strong indication that the IS charmonia {rjc and 
J/ip) survive in the deconfined phase to temperatures as high as l.GTc [lH HSl |43j. A 
detailed study of the systematic effects show, however, that the reconstruction of the 
charmonium spectral function is not reliable at high temperatures [16], in particular 
the presence of peaks corresponding to bound states cannot be reliably established. 
The only statement that can be is that the spectral function does not show significant 
changes withing errors of the calculations. Recently quarkonium spectral functions have 
been studied using potential models and lattice data for the free energy of static quark 
anti-quark pair [36]. These calculations show that all charmonia states are dissolved at 
temperatures smaller than 1.5Tc, but the Euclidean correlators do not show significant 
changes and are in fairly good agreement with available lattice data both for charmonium 
[i3l HH] and bottomonium [Ml H?]. This is due to the fact that even in absence 
of bound states quarkonium spectral functions show significant enhancement in the 
threshold region [36]. Therefore previous statements about quarkonia survival at high 
temperatures have to be revisited. The large enhancement of the quarkonium correlators 
above deconfinement in the scalar and axial-vector channel can be understood in terms 
of the zero mode contribution [361 SH] and not due to the dissolution of the IP states as 
previously thought. Similar, though smaller in magnitude, enhancement of quarkonium 
correlators due to zero mode is seen also in the vector channel [16]. Here it is related to 
heavy quark transport [ll[35]. In the vector channel the spectral function at very small 
frequency, i.e. the transport contribution is given by ^ 

c^u'"i^) = X,iT)vLr^l^J^^, (5) 

with Xq{T) being the quark number susceptibility for charm or bottom quarks and 
Vtherm IS their thermal velocity. The width of the low energy transport contribution of 
the spectral functions is rj = T/ (mgD) with nig being the heavy quark mass and D is the 
heavy quark diffusion constant [1]. Since T/rUg is large rj is quite small and the transport 
peak in the spectral functions gives an almost constant contribution to the correlator. 
Deviation from the constant gives information about the value of D. However, existing 
lattice data are not precise enough to constrain its value. Existing lattice data can 
provide information about the thermal velocity of heavy quarks. On the lattice one 
can also calculate the temporal component of the vector correlator which is stricly 
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Figure 2. The thermal velocity of heavy quarks estimated using anisotropic lattice 
data [15] (filled symbols) and isotropic lattice data 43J (open symbols). 



constant Goo(t) = —Xq(T)T |1] because of charge conservation. Therefore calculating 
the transport contribtion in the vector channel and deviding it by the temporal vector 
correlators gives an estimate of the thermal velocity of the heavy quarks. In Fig. [2] I 
show the estimated thermal velocity of the charm quarks using lattice data from isotropic 
and anisotropic lattices. 

Note that for free quarks the thermal velocity squared is just T/rUg. As one can 
see from the figure the lattice data indicate significantly smaller thermal velocity even 
assuming a fairly large value for the effective quark mass rUc = 1.7GeV. Furthermore, we 
see a dramatic decrease of the thermal velocity at temperatures close to the transition 
temperature. 

Although the above mentioned lattice studies were performed in quenched 
approximation we do not expect the picture to change when dynamical quarks are 
included in the calculations since recent calculations in 2-flavor QCD show very similar 
temperature dependence of charmonium correlators |49j . 

The spectral function for light mesons as well as the spectral function of the 
energy momentum tensor has been calculated on the lattice in quenched approximation 
[501 ISTl [52l [53] . However, unhke in the quarkonia case the systematic errors in these 
calculations are not well understood. 

4. Summary 

Significant progress has been achived in lattice calculations of thermodynamic quantities 
using improved staggered fermions. Pressure, energy density and entropy desnity 
can be reliably calculated at high temperatures when improved actions are used. 
Different lattice calculations show that for the physical quark masses the transition 
to the deconfined phase is not a true phase transition but a crossover. There is some 
controversy, however, concerning the location of the crossover. Lattice calculations 
provide detailed information about screening of static quarks which is important for 
the fate of heavy quarkonia in the quark gluon plasma. Some progress has been made 
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in calculating spectral functions on the lattice, however, much more work is needed in 
this case. One interestimg result have been obtained for the thermal velocity of heavy 
quarks. 
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